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Abstract 

We study excitations of the local field (locsitons) in nanoscale two-dimensional (2D) lattices of 
strongly interacting resonant atoms and various unusual effects associated with them. Locsitons 
in low-dimensional systems and the resulting spatial strata and more complex patterns on a scale 
of just a few atoms were predicted by us earlier [A. E. Kaplan and S. N. Volkov, Phys. Rev. Lett. 
101, 133902 (2008)]. These effects present a radical departure from the classical Lorentz-Lorenz 
theory of the local field (LF), which assumes that the LF is virtually uniform on this scale. We 
demonstrate that the strata and patterns in the 2D lattices may be described as an interference of 
plane- wave locsitons, build an analytic model for such unbounded locsitons, and derive and analyze 
dispersion relations for the locsitons in an equilateral triangular lattice. We draw useful analogies 
between one-dimensional and 2D locsitons, but also show that the 2D case enables locsitons with 
the most diverse and unusual properties. Using the nearest-neighbor approximation, we find the 
locsiton frequency band for different mutual orientations of the lattice and the incident field. We 
demonstrate a formation of distinct vector locsiton patterns consisting of multiple vortices in the 
LF distribution and suggest a way to design finite 2D lattices that exhibit such patterns at certain 
frequencies. We illustrate the role of lattice defects in supporting localized locsitons and also 
demonstrate the existence of "magic shapes", for which the LF suppression at the exact atomic 
resonance is cancelled. 

PACS numbers: 42.65.Pc, 85.50.-n 
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I. INTRODUCTION 

In our recent Letter [l| we predicted nanoscale field patterns (stratification) emerging in 
one-dimensional (ID) arrays and two-dimensional (2D) lattices of strongly interacting atoms, 
driven by a radiation nearly resonant to the atomic transition. We predicted excitation of 
so called locsitons and a host of related effects. A general formulation of the problem and a 
more detailed theory for ID arrays was presented in our most recent paper j2|. The present 
paper is an extension of toward the theory of 2D lattices of resonant atoms, which 

produce a much richer set of effects. We construct here a detailed theory of interactions 
in the system by developing different 2D versions of the nearest-neighbor approximation 
(NNA), including the "near-ring" approximation (NRA). We also derive dispersion relations 
for various lattice-polarization configurations for all locsiton wave vectors within the cor- 
responding first Brillouin zones. Our theory predicts such phenomena as subwavelength 
multicell patterns, including multivortex locsiton excitations, and locsitons localized near 
lattice defects. Further on, we predict "magic shapes" of nanosize groups of atoms, which 
reverse the effect of a resonant locsiton suppression present in all but few configurations. 
The simplest magic configuration which can be cut out of a triangular lattice is a six-point 
star with an atom at its center, which makes the lowest "magic number" of atoms to be 13. 

The predicted effects would be totally unexpected within the standard theory of local 
fields 3|] going back to the works of Lorentz jj] and Lorenz Sj]. That celebrated theory 
asserts that the microscopic electric field El acting upon any given atom in a medium — the 
local field (LF) — differs from the macroscopic field E of the electromagnetic wave, because 
electric dipoles induced in neighboring atoms produce extra field to supplement the field 
of the incident wave. This difference is significant in dense media, where the interatomic 
interactions are sufficiently strong. Under such conditions, typical interatomic distances are 
much shorter than the optical wavelength, and the dipole-dipole interactions between atoms 
can be treated as quasistatic. The major point of the Lorentz-Lorenz theory (LLT) of local 
fields is that the LLT contains a fundamental assumption (which often remains implicit 
and unspoken in the literature) that the LF varies very insignificantly between neighboring 
atoms, much like the applied optical field on the subwavelength scale. Unsurprisingly, that 
theory results in the LF being proportional to the macroscopic electric field. El = E(e+2)/3, 
where e is the dielectric constant of the medium. 
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As we have shown in , this assumption of the LF uniformity is not universally 
applicable; moreover, it completely falls apart when the uniformity of the atomic lattice 
is disturbed by impurities, boundaries, etc. Indeed, when the interatomic interactions are 
sufficiently strong and the system is not very large, highly nonuniform LF distributions 
emerge, resulting in a strong stratification of the LF and atomic excitations. This effect 
is best manifested in small-scale ordered arrays and lattices of atoms at near-resonance 
conditions, which allow to attain high interaction strength between neighboring atoms and to 
easily control it by tuning the laser frequency. We have shown [l, 2] that when the interaction 
with neighboring atoms becomes comparable to that with the external field, so that the 
interaction strength exceeds some critical value, the system will support LF excitations, 
which we call locsitons. In finite-size arrays and lattices, standing waves of locsitons may 
form nanoscale strata and complex patterns in the LF (and hence, in the induced atomic 
dipoles). A typical spacing between atoms in the arrays and lattices exhibiting the LF 
stratification is a few orders of magnitude shorter than the wavelength of light, so the 
quasistatic approximation of the standard LLT can still be used. 

It is worth noting that locsitons are basically a linear phenomenon and can be excited by 
a weak incident field. We want to stress that locsitons, i.e., spatially nonuniform solutions, 
are not new stable-state alternatives to a presumably unstable uniform Lorentz solution at 
certain interaction parameters; the stability or instability of the solution is not an issue here. 
The emergence of locsitons is determined by the boundary conditions in a finite structure, so 
a locsiton is essentially the only physical solution, which thus replaces the uniform Lorentz 
solution. Absorption plays an important role here, as it directly affects spatial attenuation 
of locsitons and thus the maximum distance to a boundary, defect, or other inhomogeneity 
where locsitons can appear. In particular, as we have shown in ^ , the size of an array that 
can support well pronounced locsitons is directly related to the characteristic absorption 
length. Locsitons vanish in the bulk of a crystal sufficiently far away from boundaries or 
defects. 

Due to recent advances in fabricating nanoscale structures, the observation and practi- 
cal applications of the LF nanostratification are becoming a reality. Theoretically, strongly 
interacting resonant particles discussed in \i] do not have to be atoms, but may also be 
quantum dots, molecules, clusters, etc. However, one has to remember that one of the ma- 
jor conditions for the structure to support the locsitons is that the interaction strength has 
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to exceed a critical value. Since this strength is proportional to the square of an individual 
dipole momentum and inversely proportional to the cube of the interparticle spacing and 
the linewidth of the particle resonance (see below), the atoms may become preferred can- 
didates. A very high finesse of atomic resonances (i.e., their narrow linewidth), compared, 
for example, to plasmons (see below), also contributes greatly to the phenomenon, allowing 
one to see high-order locsiton resonances. 

To observe locsitons, one have to create conditions to couple them efficiently to an op- 
tical or some other kind of a probe. In we suggested a few promising methods of 
locsiton detection. In particular, locsitons could be observed via size-related resonances in 
a scattering of laser radiation or via x-ray or electron-energy-loss spectroscopy. 

There are many potential applications of locsitons; here we will mention two of them 
which were discussed in [1]. It was shown in jl| that in the presence of a sufficiently strong 
optical field (i.e., in the nonlinear case), the LF in ID arrays of strongly coupled dipoles 
can exhibit optical bistability, which could be used to design nanoscale all-dielectric logic 
elements and switches. Such devices might complement currently used semiconductor-based 
electronic circuits. Another potential application of locsitons could be based on the extreme 
sensitivity of size-related locsiton resonances to the size and shape of the system. At the 
exact atomic resonance, the field is normally "pushed out" of the atomic system, unless it 
has a certain "magic shape" [l| . Consequently, such "magic structures" of atoms could find 
applications in designing nanoscale biosensors. 

Removing the assumption of the LLT that the dipoles in the medium oscillate in lock- 
step with the incident electromagnetic wave is a substantial paradigm shift in the theory of 
light-matter interaction. Locsitons predicted within our broader approach are a new phe- 
nomenon, although we can provide some incomplete but illustrative analogies from other 
areas of physics. For example, short- and long-wavelength strata in [U] are reminiscent of 
ferromagnetic and antiferromagnetic arrangements of static magnetic dipoles in the Ising 
model. The LLT is, on the other hand, more similar to the mean-field approach of the 
Curie- Weiss theory for magnetic media [6[. The Ising model is known to have richer con- 
sequences than the Curie- Weiss theory. Our case is, however, substantially different and 
most of all, more versatile than the Ising model. Indeed, instead of being static, as in the 
Ising model, the atomic dipoles are induced by the applied optical field and can oscillate 
with arbitrary amplitude and phase. By their nature, locsitons may be classified as Frenkel 
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excitons j?], because there is no charge transfer between atoms and the dipole interaction 
is due to bound electrons. Some of the locsiton effects, first of all, wave resonances, may 
be viewed as analogues of other types of oscillations and waves in condensed matter, like 

fin 

plasmons and phonons [7.], as well as low-dimensional effects, like surface plasmons [8|, |9[, 
size-related resonances in thin metal films [l^ or long organic molecules pjj], "quantum 



carpets" 



12|, arrays of pendulums or electronic circuits etc. Approaching from another 



perspective, one can view the formation of a locsiton band as a Rabi broadening of the 
atomic resonance due to strong interatomic interactions, which is essentially similar to the 
band formation in solid-state theory. However, although all waves and oscillations may be 
said to have something in common, locsitons form a distinct new class of phenomena be- 
cause of the nonconductive, dielectric, nature of their optical response and a strong coupling 
between atoms, which is needed for attaining dramatic size-related resonances, magic con- 
figurations, and other interesting effects. A separate issue outside of the scope of this paper 
is how the structures can be fabricated or arranged. One can envision placing the atoms in 
a controlled way on the surface of suitable dielectric materials; recent developments in the 
atomic- and ion-traps technology allow for arranging atoms in vacuum as ID arrays and 2D 



"crystals" in the so called wire traps 



13|. 



Our paper is structured as follows. In Sec. [TTlwe outline our problem and present some 
general formulas; for more details the reader should refer to our recent paper |2|. In Sec. IIIII 
we describe locsitons in infinite, unbounded, 2D lattices of resonant atoms in the case when 
the incident field is polarized in the lattice plane. The equations for the LF and the dispersion 
relations for the locsitons are first obtained in the NRA and then in the more precise NNA. 
In Sec. HV] we discuss effects that arise in finite 2D lattices, in particular, formation of 2D 
patterns of the LF, due to size-related locsiton resonances, and the "magic" cancellation 
of the resonant LF suppression. In Sec. |V] we describe 2D locsitons in the case when the 
incident field is normal to the lattice plane. Sec. I VI I summarizes main results of our paper. 



II. 2D LATTICES: GENERAL MODEL 



Let us consider a 2D lattice of strongly resonant identical particles, which we will further 
call "atoms". We will assume that, for the incident laser frequency u near their resonant 
frequency cjq, these atoms can be described by a two- level model with the transition dipole 
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moment da- In the linear case, i.e., when the laser intensity is significantly lower than the 
saturation intensity, the results for the two-level model coincide exactly with those for the 
classical harmonic oscillator, see We further only consider lattices of atoms interacting 
via quasistatic near-field dipole forces. In general, this assumption is valid if, on one hand, 
the minimum separation la between atoms is not too small, so that their atomic orbitals 
do not overlap, and, on the other hand, the lattice is finite and its overall dimensions are 
smaller than the laser wavelength A. However, in most cases, in particular within the NNA, 
which is the most common one and used throughout this paper, the sufficient condition is 
much less restrictive, only requiring that the interatomic separation ^ A, which is the 
same limit as in the standard LLT [3|. 

The LF acting upon a given atom located at a point r is a superposition of the incident 
("external") field Ejn and the sum of the fields Edp(r, r') from surrounding dipoles at all 
other lattice positions r' over the entire lattice: 

EL(r) = Ei,(r) + Edp(r,r'), (1) 

lattice 

where the near fields of the surrounding atomic dipoles p(r') are dominated by the non- 
radiative (quasi-static) components |14| . 

Eap(r.r-)^ ^"IP';|>;_"l^-P''-'' . (2) 

Here u = (r — r')/|r — r'| is a unit vector along the line connecting the two atoms, e is the 
background dielectric constant (e = 1 in vacuum and e 7^ 1 if a host medium is present). 
The atomic dipoles p(r') are, in turn, induced by the LF acting upon them. In the case of 
linear optical response of a two-level atom, its induced dipole moment is 

where 6 = T{u — uo) is a dimensionless laser frequency detuning and T = 2/r is the 
transverse relaxation time of the atom with a resonant homogeneous linewidth F. The 
condition that the electron orbitals of neighboring atoms do not overlap implies that la 3> 
\da\/e and justifies our use of a semiclassical approach. Nonlinear effects in systems of 
resonant atoms could be included into this picture via the saturation nonlinearity of the two- 
level system. As we have recently demonstrated [l, 2], the nonlinearity enables interesting 
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effects with promising applications, like the optical bistability and hysteresis. In the present 
paper, however, we will only consider linear effects in 2D lattices of resonant atoms. 
We can rewrite Eqs. ([I])-(I3]) in the following closed form, 

E,(r) = E,„(r) - ^ ^ 

lattice I I 

x{3ii[EL(r')-u]-EL(r')}, (4) 

where we introduced the dimensionless strength Q of the dipole-dipole coupling between 
neighboring atoms, which is easily controlled through the normalized frequency detuning S, 

The maximum normalized strength of interaction between neighboring atoms, 

is reached at zero detuning of the incident laser frequency from the atomic resonance (6 = 0). 
Note that this is the point where our approach drastically departs from the conventional LF 
theory, where the LF is implied quasi-uniform on a scale of many la, so that EL(r) ^ EL(r'). 
More details on the general LF model that we use in this paper, including some quantitative 
estimates of Qa for realistic systems, can be found in our recent pubhcation 

In the present paper we restrict our considerations to the nearest-neighbor approximation 
(NNA), in which only interactions between the closest atoms are taken into account in 



Eq. (j4]). As we demonstrated in , this approximation leads to qualitatively similar 
results compared to the full solution, that takes into account interactions between all pairs 
of atoms in the system. The NNA allows us to derive simpler analytic expressions and 
undertake numerical simulations for reasonably large lattices. 



III. LOCSITONS IN TRIANGULAR LATTICES: IN-PLANE POLARIZATION 

In this section we consider the emergence and properties of locsitons in infinite unbounded 
planar lattices of atoms. Here we aim at studying locsitons in their "simplest" form, without 
the system boundaries complicating the picture. At the same time, we will continue to 
assume the near-field character of the interatomic interaction, because our ultimate goal is 
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FIG. 1: (a) Geometry of an equilateral triangular lattice of resonant atoms, (b) The corresponding 
reciprocal lattice, shown with open circles, and some high-symmetry points and directions in the 
first Brillouin zone. The dashed lines in (b) illustrate constructing the first Brillouin zone as the 
Wigner-Seitz cell of the reciprocal lattice. 

to study systems that are much smaller than A. For the same reason, we will assume the 
incident field to be uniform throughout the lattice, Ein(r) = Ein, although our formalism is 
not limited to this case. As we have proven in [l, 2| for linear (ID) systems of interacting 
atoms, such an approach is indeed of a great help for understanding the locsitons' behavior. 

We consider an equilateral triangular lattice (also known as a hexagonal lattice) of atoms, 
which has a six-fold rotation symmetry and belongs to the plane symmetry group (wallpaper 
group) pdm. [See Fig. [I^a).] This lattice type has a remarkable property of providing the 
most close-packed configuration of identical circular objects in a plane. Each atom in the 
lattice has six neighbors at the distances of la- The "second layer" of neighbors is removed 
by y/Sla or 2/^; considering the fast decrease of the dipole-dipole interaction force with the 
distance (oc l/r^), the NNA which ignores interaction with the second and further layers of 
neighbors is expected to work well. 

n n 

As in the case of ID arrays of atoms [1, 2], two major cases can be studied separately, 
where the incident laser field Ein hes in plane (the "||" case) or is normal to the plane of 
the lattice (the "_L" case). The linear optical response in the case of any other incident 
polarization can be obtained using a superposition of these two configurations. Of the 
two major cases, the "||" configuration is by far more interesting, exhibiting richer locsiton 
behavior and differing significantly from the ID case due to the crucial in-plane anisotropy 
of the dipole-dipole interaction. It is also easier to implement, as the laser beam may be 
incident at the normal to the lattice of atoms, with the incident field being virtually uniform 
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on a scale of many wavelengths. The "_L" configuration presents a much closer analogy to 
the ID problems considered in , although some details inevitably differ. We will discuss 
it below in Sec. El Further in this section, we only consider the "||" case, where both Ejn 
and El lie in the lattice plain. 



A. Near-ring approximation 

Optical response of a triangular lattice of atoms is, in general, anisotropic: it depends 
on the orientation of Ein in the lattice plane with respect to the symmetry directions of 
the lattice. We have demonstrated, however, that one can use the near-ririg approximation 
(NRA) to describe the behavior of long- wavelength locsitons in the lattice jl| in an isotropic 
fashion. In this approximation, the contribution from the surrounding dipoles to the LF El 
[the second term in Eq. ([1])] is substituted with a field from an effective dipole ring with the 
radius /„, such that the polarizability of the nearest six atoms is evenly redistributed along 
the ring. The NRA is thus a further simplification of the NNA, making the model isotropic 
in the lattice plane. Within the NRA we then replace the summation in Eq. (jl]) with an 
integration over the imaginary ring, so that the equation for the LF becomes 

EL(r) = Ei, - r EL(r + /,u) (3 cos^ 9 - 1) d9, (7) 

where u is a unit vector in the direction from the center r of the effective ring to a point at 
the ring, and 9 is the polar angle of u counted from the direction of Ein in the lattice plane. 
The strength of the dipole-dipole coupling between neighboring atoms is still given by the 
dimensionless parameter Q [Eq. ([5])]. 

Equation (JTj) has a uniform (Lorentz) solution 

= T — T^— — Ein, (8) 

^ - Oll + « 



where 



Sl = -iQa (9) 



is the normalized frequency detuning at which the Lorentz-Lorenz resonance is achieved in 
the triangular lattice of atoms within the NRA. As we will see in Sec. El 5ll is different 
for the "II" and "_L" configurations, but to simplify the formulas, we omit the index "||" 
everywhere in this section and in Sec. ITVl except for Eq. (]9]). From Eq. (|8]) one can see that, 



if ~ Qa ^ 1, the uniform Lorentz LF is suppressed when the laser is tuned closely to 
the exact atomic resonance, 5 ~ 0, reaching its minimum intensity 



lEi, 



|2 



The LF in this case is effectively "pushed out" by the lattice atoms. Interestingly, a huge 
LF enhancement is reached at a red-shifted frequency, at 5 ~ 5ll < 0, where 

|EL|Lx^(l + '5LL)|Einp. (11) 

Note that Eqs. ([8]), (fTOj) . and (fTT!) are very similar to the corresponding equations for 
ID arrays of atoms which were discussed in [l|, |2[. The only difference is that the relation 
between and Qa there was 5ll = ~Qa if ^in is parallel to the array (and dipoles are 
aligned "head-to-tail") and 5ll = Qa/2 if Ejn is perpendicular to the array (and dipoles are 
aligned "side-to-side"), assuming the NNA [see, e.g., Eq. (3.4) of Ref. |2l]. Here, within 
the NRA, 5ll does not depend on the incident field polarization in the lattice plane. Quite 
naturally, its value (|9]) lies in-between the two values for the ID array, because mutual 
orientations of different pairs of dipoles in the lattice vary between the two extremes. 

The frequency dependence of a spatially uniform LF, like in Eq. ([H]), and the associated 
Lorentz shift are long known phenomena. Similar effect was also observed experimentally 



in alkali vapors [15|. The unusual new phenomenon is that in ordered low- dimensional 
structures there are spatially varying solutions, which we call locsitons [1], that emerge at 
some values of Q in addition to the uniform LF. We will look for the locsitons in the form 
of 2D plane-wave excitations of the LF: 

AELOcexp(iq-r/g, (12) 

where q is the normalized wave vector of the locsiton. By substituting El = El + AEl into 
Eq. ([7]) we obtain the dispersion relation for the wave vector q in an integral form: 

1 + — r(3 cos 29 + 1) cos[g cos(^ - i;)] dO = 0, (13) 

4:71 Jo 

where ip is the polar angle of q counted from the direction of Ei^ in the lattice plane. Using 
the standard expansion of trigonometric functions with a harmonic argument into Bessel 
functions [see, e.g., Eq. (21.8-25a) in Ref. [igI], 

oo 

cos{qsm(j)) = Jo{q) + 2 ^ J2m(g) cos(2m0), (14) 



m=l 
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where Jm{<i) is the Bessel function of the first kind, we evaluate the integral in Eq. (fT3!) in 
an explicit form and write Eq. f lT3|) as 



By substituting Q from Eqs. ([5]) and ([9]), we can rewrite the last formula as an explicit 
dispersion relation connecting the normalized detuning 5 with the normalized locsiton wave 
vector q which is represented by its polar coordinates q and ip: 



In the limit of low absorption, i.e., ^ 1? the r. h. s. of Eq. f lT6|) may be replaced with 
a more detailed analysis of the dispersion relation shows that spatial oscillations — 
locsitons with almost real q — emerge in a limited frequency band around ujq with a normal- 
ized bandwidth ~ I^llI- It is exactly the frequency range where the interaction between 
resonant atomic dipoles may become comparable to or much stronger than the effects of 
the external field Ejn. The dispersion dependence = -D^^^(q) in this case is shown 

in Fig. Et^a), where the external field is assumed to be aligned with the x axis. At every 
given laser frequency w, locsitons with a whole range of wave vectors q may be excited. 
This set of q, all having different orientations, is represented by an isoline for the respective 
^/^■LL = -D^^'^(q) in Fig. [2]^a). Recalling that 5ll < [see Eq. Q], we may notice that 
solid isolines (red shaded areas) correspond to red-shifted laser frequencies {u < uq), while 
dashed isolines (blue shaded areas) correspond to blue-shifted frequencies. The non-circular 
shape of the isolines reflects the highly anisotropic dispersion dependence for locsitons in 
the 2D lattice in the case of short-wave locsitons. In particular, at any given normalized 
frequency detuning 6, locsitons with different orientations of q with respect to Ein may have 
very different q and, consequently, different wavelengths. At the same time, because the 
NRA describes the lattice in an averaged way, the orientation of Ejn with respect to the 
lattice does not affect the result. 

In the general case, especially when looking at the limitations on the size of the structure, 
one needs to consider complex locsiton wave vectors q = q' + iq" in the dispersion relations 
(fT3ll or (fT6l) . like it was done for ID arrays of resonant atoms in This would allow to 
describe the dissipation of 2D locsitons, which is most prominent near the Lorentz resonance, 
and evanescent locsitons, which exist outside the locsiton band. One can also calculate the 




(15) 




'LL 



(16) 
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FIG. 2: (Color online) Dispersion dependences for locsitons in a triangular lattice of atoms (a) 
in the near-ring approximation; (b) for Ein || TK; (c) for Ein || TM. Isolines of equal -D2(<l) are 
spaced at 0.25 (i. e., they correspond to the ticks on the color bars); solid isolines and red shading 
correspond to positive D2(q)) dashed isolines and blue shading correspond to negative ones, while 
the thicker solid isoline marks the zero level. The hexagons in plots (b,c) show the boundaries of 
the first Brillouin zone. 
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group velocity of the locsiton, Vgr = {la/T){dS/dq), by taking a derivative of Eq. (fT6|) . 
Due to the anisotropy of Eq. f|T6l) . one may expect that the resulting dissipation and group 
velocity are highly dependent on the orientation of q with respect to Ejn. 

By its nature, the NRA only works well for long- wave locsitons, i. e., those with relatively 
small q. Indeed, for a long-wave locsiton, the LF at neighboring atoms differs insignificantly, 
so one may reasonably expect that replacing the neighboring dipoles at their actual positions 
with an effective dipole ring will not cause a significant change to the result. Our results 
obtained without resorting to the NRA (see Sec. IIIIBI) . suggest that the NRA provides a 
good quantitative description of locsitons for q ^ 7r/2, while it gives reasonable qualitative 
estimates for g up to ~ tt. One may use the two first terms from the Taylor series expansion 
to analyze Eq. f lTB]) at small q, which is easy to d o hj recalling that Jo{q) = 1 — -|- 0{q^) 



and J2(g) = |r + O(g^) [cf. Eq. (21.8-3) in Ref. [I6l], so that 



D^RA^q) ^ 1 - -g'[2 + 3 cos(2^)] + O(g^). (17) 
8 

In the long- wavelength limit, D^^^(q = 0) = 1 is reached at S = ^llj yielding the uniform, 
Lorentz, solution for the LF. The second term in the r. h. s. of Eq. f lT7|) indicates that the 
anisotropy with respect to the locsiton polarization shows up even at very small q, where 
Eqs. f lT6|) and f lT7|) can be rewritten as 

2 + 3cos(2V') 5ll ■ 

Note that in describing locsitons at very small q (i. e., for 6 ~ 5ll), the imaginary unit in 
the last formula cannot always be neglected compared to 6 — The locsiton dissipation 
in this case may become significant, depending on and the necessity of considering 
complex q could appreciably complicate calculations. We will not go into the details of this 
special case in the present paper. 

In ID arrays of atoms, the Lorentz resonance, 6 = always coincides with one of the 
edges of the locsiton frequency band In that case, when one tunes the frequency 

of the incident laser beam towards Uq from the side of the Lorentz resonance, locsitons 
with the longest wavelength emerge first very close to the Lorentz resonance. The locsiton 
wavelength subsequently decreases as we approach and move past uq. In a 2D lattice, the 
Lorentz resonance lies within the locsiton band, which means that, when similarly tuning 
the laser frequency towards Uq, short-wave locsitons will be excited first. In particular, 
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locsitons with ip ~ ±7r/2 (i.e., nearly transverse locsitons) may be exited at > 1, 

and can be thus viewed as the "easiest to excite" on the Lorentz side of the band. On the 
opposite, "anti-Lorentz" side of the band, where ~ —1, nearly longitudinal locsitons 

with ^/^ ^ or TT lie closer to the band edge and thus are easier to excite. The exact positions 
of the edges of the locsiton band cannot be found within the NRA, because minima and 
maxima of D^^^(q) are reached at such q where the NRA may only be used for qualitative 
estimates. The case of larger q is addressed when we go beyond the NRA in the next 
subsection. 

B. Locsitons in the first Brillouin zone 

When going beyond the NRA, the orientation of Ejn within the lattice plane becomes an 
important factor, except for small q. Staying within the NNA, i. e., only taking into account 
the six nearest neighbors in Eq. (jl]) (but individually, instead of them being washed out 
over the ring, as in the NRA), we are still able to approach the problem analytically. The 
resulting equation is 

EL(r) = Ei,-^^{3uK[EL(r + /„UK)-UK] 

UK 

-EL(r + /,UK)}, (19) 

where Uk denotes any of the six unit vectors pointing in the directions from the atom (located 
at r) to one of its nearest neighbors. 

The uniform, Lorentz, solution of Eq. ( [T9|) is still given by Eq. (|8]) and ([9]), which supports 
the above-mentioned convergence of the NRA and NNA results at g — 0. Spatially varying 
locsiton solutions are found as in the previous subsection by using the ansatz (JT2l) in Eq. (fTOll . 
The corresponding dispersion relation 6{q) for locsitons in a 2D triangular lattice can be 
now written as 

^^'^^^(q) ^ ^ f cos 26^ + ^) cos[g cos(0„ -^)] = thl_ (20) 

where 6n = Oo + nir/S, and q is represented by its polar coordinates q and ip. The orientation 
of the lattice with respect to the incident field Ein is described by 6*0, which is the angle 
that one of the vectors uk makes with Ejn (the fact that 6q is not unique does not affect the 
result). The ultimate proof that the results of the NRA and the more precise NNA converge 
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for long locsiton wavelengths can be obtained by taking a Taylor series expansion of Eq. ( 120|) 
at g — 7- 0. By only retaining terms up to the order of and calculating all the necessary 
sums and products of trigonometric functions, we find that Eqs. (fT7|l and ( ITS]) still hold in 
the NNA. 

It is sufficient to find locsitons with q lying within the first Brillouin zone of the reciprocal 
lattice, because any solutions with q lying outside the first Brillouin zone are physically 
equivalent to them due to the discrete nature of our system. The dispersion relation D^^'^(q) 
in Eq. ( 120|) also has the required symmetry and periodicity. It is common to use high- 
symmetry points in the first Brillouin zone to denote most interesting directions in the 
lattice [see Fig. [U^b)]. Note that in terms of uk we may write FK || uk and FM ± uk. 

The dispersion dependence = -D^'^^(q) in the case of ^ 1 is shown in 

Fig. |2]|^b,c) for two different orientations of Ejn with respect to the lattice: in Fig. [2]|^b) 
^0 = (i.e., Ein II FK), while in Fig.EJ^c) 9q = tx/2 (i.e., Ei^ || FM). The x axis on the plots 
is aligned with the external field Ein; the hexagonal boundaries of the first Brillouin zones 
are shown with thicker dashed lines. The central parts of all three plots in Fig. [2] are very 
similar, which refiects our finding that Eqs. ( 1T7|1 and ( fT8|) hold for long locsiton wavelengths 
within both NRA and NNA. At the same time, significant differences accumulate closer to 
the boundaries of the first Brillouin zone. Fig. [2]^b,c) also show that for both orientations of 
Ein the maxima and minima of D2^^{^q) are attained at the zone boundaries, specifically, 
at different M points. 

To facilitate the comparison of the three plots in Fig.|2]and finding the edges of the locsiton 
band, the NNA dispersion dependencies for two high-symmetry directions are presented in 
Fig. |3] and compared to the the NRA result. For both longitudinal and transverse locsitons. 
Fig. [3] shows corresponding cross-sections of the three plots of Fig. |2j As we already noted 
earlier in this paper, the NRA result is only meaningful for g up to ~ vr; Fig. [3] suggests 
that it is a good approximation for the NNA result at g ^ 7r/2, regardless of the orientation 
of Ein with respect to the lattice. It is instructive to give explicit analytic expressions for 
each of the four NNA-based dependencies in Fig. [3l The respective dispersion relations 
derived from Eq. fl20|) and the corresponding ranges for obtained in the assumption 

that |(5ll| ^ 1, are as follows: 
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FIG. 3: Dispersion dependences for (a) transverse locsitons (q _L Ein) and (b) longitudinal loc- 
sitons (q || Ein) in a triangular lattice of atoms. Different curves correspond to the near-ring 
approximations (NRA) and the NNA for two different orientations of q with respect to the lattice. 



(a) Ein II TK, q II TM ± Ei^ [long-dashed line in Fig. ^a)] 

S_ 

5ll 



COS — — = 4 — 3— — , [21) 



5ll 3 

(b) Ein II FM, q II FK ± Ejn [solid line in Fig. El^a)], 



1 < < 4; (22) 



cos I 4 (5 ± ^57- 48^ j, (23) 

-4^A^^4' (24) 
3 Oll Id 

where the minimum of is reached outside the first Brillouin zone; this point more 
appropriately belongs to the case (d) below; 
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TK, q II TK II Ein [solid line in Fig. Mh)] 
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where the maximum of is reached outside the first Brillouin zone, which point 

more appropriately belongs to the case (a); 



(d) Ein II TM, q II TM || Ei^ [long-dashed line in Fig. Mh)], 



cos 



qV3 1 



5^ + ^ 



6 + i 



1 6 

2- < 

3 - 5ll 



< 1. 



(27) 



(28) 



The inequalities fl21|) and also represent the edges of the locsiton bands for the cases 
of Ein II TM and E^ || FK, respectively (which explains our desire not to restrict the range 
of q to the first Brillouin zone when obtaining these inequalities). Note that locsitons with 
the largest red shift (5/5ll > 0) can be achieved with Ein || TK, while locsitons with the 
largest blue shift (5/5ll < 0) can be achieved with Ein || TM. Therefore, re-orienting the 
lattice with respect to the polarization of the incident laser beam may assist in controlling 
the type of locsitons excited in the lattice. 



IV. FINITE LATTICES: IN-PLANE POLARIZATION 

In this section, like in Sec. IIIIl we only consider the case where Ein is spatially uniform 
and lies in the lattice plane (the "||" case), which is easier to achieve if the laser beam 
is incident normally to the lattice plane. The presence of boundaries and defects in 2D 
lattices of resonant atoms can cause various locsitonic effects, including giant LF resonances. 



nese 



a, 



formation of dipole strata, and "magic" cancellation of the resonant LF suppression, 
effects are similar in their nature to their counterparts in ID arrays of resonant atoms (l 
but their manifestations are much more diverse because of the inherent anisotropy of the 
dipole-dipole interaction in the 2D case, especially in the "||" geometry. Because of this, an 
all-encompassing study of finite 2D lattices is hardly possible within this pilot study on the 
subject, so we will restrict ourselves to providing some of the most characteristic results, 
which emphasize distinctions from the ID problem. 
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A. Size-related resonances and local-field patterns 



Locsitons in a finite ID array of atoms exfiibit size-related resonances, cliaracterized by 
large increases in their amplitudes at certain frequencies within the locsitonic band, because 
locsitons are reflected at the boundaries and form standing waves (strata) [sj. Essentially, 
these resonances correspond to locsiton eigenmodes defined by the boundaries. For the long- 
wave strata they are similar to oscillations of a quantum particle in a box, as, e.g., for ID- 
confined electrons 10|, lUl , or a common violin string. It is natural to expect such resonances 



and eigenmodes to also exist in higher dimensions, in particular, in finite 2D lattices, where 
we also encounter locsiton reflections at the boundaries. An important distinction of the 
2D case is that the wave vector q of a locsiton may have an arbitrary orientation in the 
lattice plane with respect to the incident field Ei^. Multiple reflections and interference of 
locsitons with all possible q quickly make the whole picture very complicated and highly 
susceptible to minor changes to the size and shape of the lattice patch. We found that 
at certain geometries only a limited number of locsiton eigenmodes are dominant. Their 
interference produces various dipole patterns and strata; some of them are reminiscent to 
"quantum carpets" An important issue is, therefore, how one can control the locsiton 
patterns via the geometry of the lattice patch and the frequency and polarization of the 
laser beam. 

One way to engineer a distinct 2D locsiton pattern is to start with a rectangular lattice 
patch and ensure that size-related resonances are achieved for locsitons with wave vectors 
parallel to its boundaries. We have to choose the lattice shape, such that the size-related 
locsiton resonances emerge in both dimensions at the same frequency detuning 6. To simplify 
our task, we will consider long-wavelength locsitons, which are not too sensitive to the system 
sizes and thus are easier to control, and, incidentally, also form more pronounced patterns 
and are described by the simpler formula f|T8|) . 

In the limit of long- wavelength locsitons (g <^ 1, 5 ~ the dispersion relations in the 
cases (a) and (b) described in Sec. IIIIBI coincide with each other: 

= 8 ^^-^^^ + ^ (29) 

[cf. Eq. ([TSjl aX = n / 2] . In a similar manner, one obtains approximate solutions for the 
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FIG. 4: Vortices in the distribution of the local field in a nearly square patch of a 2D triangular 
lattice of atoms at (5 = —1000 and 5ll = —987.375. To avoid overcrowding of the plot, only one 
of each nine dipoles is shown. The incident light wave is polarized in the lattice plane along the 
diagonal of the lattice patch, its field is shown with a big arrow. 

cases (c) and (d), for which = 



By combining the cases (a) and (b) or the cases (c) and (d), we can achieve simultaneous 
size-related resonances, i.e., excitations of eigenmodes in both orthogonal directions in a 
patch of the 2D triangular lattice at the same laser frequency, if the patch is approximately 
square in shape. Resonances of the same order are hereby attained for locsitons with wave 
vectors pointing in the two orthogonal directions; a sufficient "squareness" of the lattice 
patch can be achieved by choosing its size (i. e., the numbers of atoms in the two directions). 
Locsitons with shorter wavelengths or with wave vectors pointing in different directions will 
be also present, but they will have no significant infiuence on the emerging dipole pattern 
due to their nonresonant nature. 

Fig. m shows vector patterns that are formed by the atomic dipoles induced by the LF. 
The atoms are arranged in a 48 x 56 patch of an equilateral triangular lattice, which results 
in approximately equal sides of the patch. The field of the incident electromagnetic wave 
is uniform and polarized along the diagonal of the patch and its frequency is close to the 




(30) 
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electronic resonance of the two-level atom. The frequency of the incident wave is so chosen 
that the third size-related resonance (in the order of increasing wavenumbers, counting only 
those resonances allowed by the symmetry of the problem; for more detail see j2|) is excited 
in each dimension. Eight distinct vortices of the LF are visible in the plot. Fig. IHonly shows 
the imaginary parts of the complex field amplitudes, because they are dominant for each of 
the resonant locsitons. (We would like to note that a pair of vortices, apparently consistent 
with 2D-locsiton patterns originated by the 1-st locsiton resonance in our classification, was 
very recently observed in numerical simulations of plasmonic excitations in a 2D lattice of 
small metallic particles jirjl.) 



B. "Magic shapes" 

As we noted in Sec. Illl Al the LF is "pushed out" of the lattice of strongly interacting 
dipoles at the exact atomic resonance [see Eq. ( |TOl) ]. This effect, which we call the resonant 
LF suppression, represents a typical LF behavior at the atomic resonance; it is not limited 
to 2D lattices, but also occurs in ID arrays of interacting atoms and in many finite 

2D structures. We have shown earlier that in the ID case, if a linear array of atoms is of a 
certain "magic size", one encounters a cancellation of the resonant LF suppression, where 
one of the size- related locsitonic resonances partially restores the LF in the system |2|. 

Finite 2D lattices and similar small systems of resonant atoms provide especially interest- 
ing examples of cancellation of the resonant LF suppression. Unlike in ID arrays of atoms, 
the "restoration" of the LF in such systems at 5 = 0, compared to that in the uniform. 



the 2D 



Lorentz, case, can be more complete (up to 100%). Like in the ID case 
"magic shapes" have a certain "cabbalistic" streak. For example, in the NNA the effect is 
most pronounced only in a system of = 13 atoms arranged as an equilateral six-point 
star with an atom at the center, for which the maximum restoration of the LF is reached, 
Ema.x/Ein ~ 1.02. The directions and relative amplitudes of the LF at the atoms in this 
system are shown in Fig. |5]^a) for Ejn || uk and in Fig. [5](b) Ein J- uk. It is very notable 
that the system is "magic" for both orientations of the incident field. One can see from the 
picture that the LF is concentrated on the outermost atoms and the one at the center, while 
the LF at the inner hexagon of atoms is almost completely suppressed. This suppression is 
a manifestation of a special case of a locsiton standing eigenwave in a finite discrete atomic 
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FIG. 5: "Magic" planar system of 13 resonant atoms for two different orientations of Ejn shown in 
plots (a) and (b). The inset illustrates the geometry of the system. 

structure. (Within the NNA, the LF zeroes out at the "empty" atomic locations outside the 
outer hexagon, so locsitons make a 2D standing wave.) In general, the zeroes (nodes) of that 
wave are located somewhere in between atoms, so at each individual atom we have a nonzero 
LF amplitude. However, in the magic atomic configuration at the precise atomic resonance, 
(5 = 0, these nodes nearly coincide with the locations of the inner-hexagon atoms. Thus, we 
have a nearly ideal picture of a 2D standing wave, with large LF intensities at the antinodcs 
(maxima), located at the central atom and the outer hexagon atoms, on one hand, and nodes 
(zeroes), located at the inner- hexagon atoms, on the other hand. To an extent, this situation 
is reminiscent of a 2D standing wave on a water surface in a round bucket with the first 
antinode at the center of the bucket, where some middle observation points are located at 
the nodes of the wave. Any symmetry distortion in this system (e. g., by attaching a foreign 
atom or molecule to it) would break the balance of the local fields in the system and bring 
back the resonant LF suppression, which is cancelled in the symmetric "magic system" . This 
effect could potentially lead to designing nanometer-scale sensors for detecting various bio- 
logical molecules, etc. For example, such a nanodevice may include target-specific receptor 
molecules that form a locsiton-supporting "magic" system. A localized locsiton then would 
get suppressed whenever a target biomolecule attaches to a receptor, otherwise the locsiton 
suppression at the electronic resonance would be "magically" cancelled. 
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V. LOCSITONS IN TRIANGULAR LATTICES: NORMAL POLARIZATION 



The "_L" configuration, where the incident field is polarized at the normal to the lattice 
plane, may be realized, for example, by creating a standing wave by two counter-propagating 
laser beams with the beam axes lying in the lattice plane. For systems much smaller than 
the laser wavelength, the incident field may be then assumed nearly uniform. Locsitons 
emerging in the "_L" configuration are similar in many respects to ID locsitons discussed in 
detail in [21], so here we only provide a brief overview of the "_L" case and outline the most 
distinctive features appearing in this geometry. 

Let us start again with describing locsitons in unbounded lattices. The most general 
equation for the LF in the "_L" case is obtained from Eq. by setting Ein(r) || EL(r) ± u: 

EL(r)=E,n(r) + | E T^E,(r'). (31) 

lattice I I 

As in the "||" case, we will assume here that the incident optical field is uniform, Ein(r) = Ein, 
which will help us to build a clearer understanding of the locsiton behavior in 2D lattices. 
For an equilateral triangular lattice, Eq. (|3T|) can be simplified using the NNA as 

ElW =Ei, + ^5]EL(r + LuK). (32) 

UK 

Equation (132|) has a uniform, Lorentz, solution, which is given by Eq. ([8]) with 5ll = ^ll 
where 

'^I^L = iQa, (33) 

which is three times the NNA value for 5ll in a ID array of atoms if Ein is perpendicular 
to the array (and the dipoles are aligned "side-to-side") j^. This is a consequence of each 
atom having now 6 instead of 2 neighbors. 

To obtain the corresponding dispersion relation 5(q) for "±" locsitons in the lattice, we 
substitute the LF EL(r) in Eq. (132|) as a sum of the Lorentz solution given by Eqs. ([8]) 
and fl55]) and 2D plane-wave excitations f[T21) . This dispersion relation can be written as 

D2^\^) = ^ E cos[g cos(^^„ - V;)] = (34) 

n=0 "^LL 

where, as in Sec. IIII B\ 0^ = Oq + rni /?>. In the long- wavelength limit, 

D,\^^(q)^l-lg2 + 0(g^), (35) 
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so that the uniform, Lorentz, sohition for the LF is reached at S = where Z?^j^^(q = 
0) = 1. The second term in the r. h. s. of Eq. f l35|) is independent of the orientation of 
q, which means that no anisotropy caused by the lattice structure is present in the long- 
wavelength limit. We may thus conclude that, compared to the "||" case, the locsitons 
in the "_L" configuration are more reminiscent of the locsitons in ID arrays of resonant 
atoms considered in There is still no complete analogy here, as, e. g., the second 

term in the r. h. s. of Eq. (I35l) differs by a factor of 1/2 from the ID result 2|. Moreover, 
dispersion relation does become anisotropic for larger q, closer to the boundaries of the 
first Brillouin zone. This anisotropy, however, is by far less pronounced than that in the "||" 
case. 

It is instructive to also obtain the dispersion relation in the NRA. By replacing the 
summation in Eq. flM|) with an integration over the "near ring", following the procedure 
outlined in Sec. IIII Al we get 

1 _ ^ r cos[g cos(^ - ^)] dO = 0. (36) 
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The resulting dispersion relation turns out to be independent of the orientation of q: 



D^l\^) ^ Mq) = (37) 



which is not surprising given the NRA applicability in the long- wavelength limit. 

While it might be somewhat harder to create a uniform incident field polarized normally 
to a 2D lattice, the resulting locsitons could be much easier to control because of the small 
anisotropy of the interatomic interactions in the "_L" geometry, compared to the "||" ge- 
ometry. For example, defects in a 2D lattice can support localized locsitons, not unlike 
the evanescent ID locsitons discussed in Compared to the complex locsiton patterns 
emerging in the "||" geometry [cf. Fig. H] these localized locsitons are more likely to form 
well-organized strata-like patterns in the "_L" geometry. 

Fig. [6] shows concentric dipole strata that are formed around a circular hole made by 
removing a few tens of atoms from a triangular lattice. The locsiton "attached" to the defect 
"decays" as the distance to the hole boundary increases, which is mostly a "diffraction" 
effect, although some contribution from the imaginary part of q (like in evanescent ID 
locsitons) is also present. In performing the numerical simulations for the plot, we made 
sure that the locsitons attached to the outside boundaries of the lattice patch (lying far 
outside the plotted region) do not interfere with the locsiton locahzed at the defect. 
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FIG. 6: Spatial strata of the normalized local field around a 15-point-wide hole in a 2D triangular 
lattice at S = fOO and 5ll = 103.5. The direction of the incident field is shown with a big arrow. 

VI. CONCLUSIONS 



In this paper we presented a detailed study of locsitons in nanoscale 2D lattices of resonant 
atoms with strong dipole interaction. These locsitons (i. e., excitations of the local field) and 
various associated effects were originally predicted in our recent publication ij. Here, we 
have built analytic models for locsitons in infinite 2D triangular lattices of atoms, based 
either on the nearest-neighbor approximation or on the simpler near-ring approximation. 

We have shown that the "in-plane" polarization geometry, where the incident laser field 
lies in the lattice plane, enables locsitons with the most unusual and diverse properties, as 

n 

compared to the ID case described in detail in [2|. In particular, the dispersion relations 
for the locsitons with an in-plane polarization are highly anisotropic with respect to the 
orientation of the locsiton polarization relative to its wave vector in the lattice plane, because 
of the highly anisotropic nature of the dipole-dipole interaction. We further demonstrated a 
method to design a finite 2D lattice, such that distinct vector locsiton patterns are formed 
at a certain laser frequency, the patterns containing multiple vortices in the local field 
distribution. 

We have also considered a remarkable effect of a cancellation of the resonant local field 
suppression jl|, which consists in the local field being able to penetrate certain "magic 
shapes" made of resonant atoms, despite the nearly-universal tendency of the local field to 
be "pushed out" of the lattice at the exact atomic resonance. In particular, we provided 
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more detail on the local field distribution in the simplest "magic shape" that can be cut out 
of a triangular lattice — a six-point star with an atom at the center. 

In the case where the incident field is polarized normally to the lattice, we found that 
locsitons bear more analogy to locsitons in ID arrays of atoms, compared to the case of an 
in-plane polarization. Finally, we illustrated the role of lattice defects in supporting localized 
locsitons. 

While this paper does not elaborate on nonlinear effects involving 2D locsitons, to be 
addressed in our future publications, we note that, similarly to the ID case O], our 
numerical simulations have shown optical bistability and hysteresis, which may be especially 
important for potential applications of 2D locsitons in designing all-dielectric nanoscale logic 
elements, devices for signal processing, etc. 
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